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Abstract 


Atmospheric structure derived from satellite, multi-channel 
radiance data is used to calculate zonal ly-averaged i^ertical motions 
in the wintertime stratosphere of both hemispheres. The Northern 
Hemisphere calculations confirm the two-celled meridional circulation 
calculated by previous authors. The Southern Hemisphere results show a 
three-celled structure with descending motion over the South Pole. The 
variability of the mean meridional circulation in the Southern Hemi- 
sphere in relation to the presence or absence of a minor midwinter 
warming is also discussed. 
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1. Introduction 


A number of studies (Reed, Wolfe and Nishimoto, 1963; Mahlman, 1966 
Julian and Libitzke, 1965) have shovm that the mean meridional circula- 
tion in the lower half of the stratosphere in the Northern Hemisphere 
winter has a two-celled structure w.' :h rising motion in low latitudes 
and over the pole and sinking motion in mid- latitudes. This zonal ly- 
averaged vertical motion pattern exists before, during and after mid- 
winter breakdowns of the polar vortex (Mahlman, 1966). In the Southern 
Hemisphere the mean meridional circulation is practically unknown. 

Rubin and Weyant (1963) do indicate a sinking motion over the pole below 
75 mb. Also Reiter (1969, 1971) on the basis of ozone measurements over 
Antarctica, postulates a general sinking motion over the role in the 
stratosphere. No evidence is available in the literature regarding the 
mean meridional circulation outside the southern polar regions. 

The primary reason for the lack of information concerning the mean 
meridional motion in the Southern Hemisphere stratosphere is, of course, 
the lack of a substantial radiosonde network in the Southern Hemisphere. 
The lack of data precludes daily analysis of the stratosphere from 
conventional data alone. In the past few years, however, data has 
become available from satellite, multi-channel infrared radiometers such 
as the Satellite Infrared Spectrometer (SIRS). Techniques have been 
developed to obtain atmospheric temperature profiles from such data (see 
for example an aritcle by Smith, Woolf and Jacob, 1970) . Using such 
temperature retrieval techniques one can obtain daily hemispheric 
analyses of temperature structure at various levels. In this study 
temperature and geopotential height information is obtained from Nimbus 
3 SIRS data. The atmospheric structure so obtained is used as a basis 
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for a computation of the zonal ly-averaged vertical motion patterns during 
winter months in both hemispheres. In the iNorthem Hemisphere the 
calculations are also carried out using the data fields of the National 
Meteorological Center (NMC). 

2. Method used to obtain atmospheric structure from Nimbus 3 SIRS 
radiances 

In this study thickness or temperature information is determined 
from the SIRS data by a regression technique. A linear, step-wise, 
least squared error, multiple regression technique is used. The depen- 
dent variables in the regression procedure are the thicknesses for the 
following layers; 1000-700 mb, 700-500 mb, 500-300 mb, 300-200 mb, 
200-100 mb, 100-50 mb, 50-30 mb and 30-10 mb. The independent variables 
are the radiances of the eight SIRS channels. Separate regression 
equations are determined for each layer and for each of the following 
latitude zones: 20-40®, 40-60® and 60-80®. 

The regression technique is based on a comparison of thickness 
information and cloud-free radiance data. When the regression -equations 
derived from the comparison are applied to other radiance data, only 
cloud-free data are used. To eliminate cloud-contaminated data from the 
comparison data set and from the application of the regression equations, 
a simple, objective "cloud check" procedure is used. The "cloud check" 
procedure is based solely on the satellite radiance information. 

The SIRS channels with weighting function peaks in the low tropo- 
sphere (channels 1 through 3) are most affected by the presence of 
clouds. Channels with weighting function peaks in the stratosphere 
(channels 7 and 8) are only rarely affected by the presence of 
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tropospheric clouds. The radiance in the window channel (channel 1, 
k=899cm“^) in the absence of clouds is related to the surface tempera- 
ture. In the presence of an overcast, the channel 1 radiance is a 
function of cloud-top temperature, which is, of course, lower than the 
surface temperature. Therefore, channel 1 radiances much below normal 
would indicate the presence of clouds. However, the radiance in channel 
1 is highly variable even in the absence of clouds because of its 
dependence on the surface temperature. Large changes in channel 1 
radiance occur along satellite tracks in the presence of sharp changes 
in surface characteristics. This is especially true along land-sea 
boundaries. Also, large diurnal changes are present in the channel 1 
radiance because of large changes in surface temperature. Channels 2 
and 3 are also affected by clouds, but are not so severely affected by 
surface characteristics as channel 1. Therefore, channels 2 and 3 are 
used for the "cloud check." 

SIRS radiance data is eliminated as being cloud-contaminated when 
the observed radiance in both channels 2 and 3 is below critical values. 
The critical values for both channels are determined by a comparison of 
SIRS data and satellite video data and are a function of latitude and 
season. Any data with both channel 2 and 3 radiances below the critical 
levels are eliminated as cloud-contaminated. The "cloud check" proce- 
dure is simple and certainly not foolproof. Some slightly cloud- 
contaminated data may still find their way into the final data set. The 
procedure as outlined above removes, on the average, about 15 percent of 
the original data points. 

To obtain the regression coefficients, comparison sets of radiance 
and thickness data are developed. The observations are near-simultaneous 
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and at, or close to, the saioe geographic position. Somewhat different 
techniques are used for the Northern Hemisphere and the Southern 
Hemisphere. In the Northern Hemisphere, for the layer 100 to 1000 mb, 
the SIRS radiance data is matched with the National Meteorological 
Center (NMC) Northern Hemisphere fields. For satellite data occurring 
within three hours of NMC map time, the NMC grid point data is linearly 
interpolated to the satellite track position. For the layers above 
100 mb, station data is matched with the satellite radiance information. 
The station location must be within 180 nautical miles of the satellite 
position, and the satellite observation time must be within three hours 
of the station observation time. 

In the Southern Hemisphere, a different technique is used. Station 
daca is matched with satellite radiance data for all levels. Because 
of the relative lack of stations in the Southern Hemisphere and the need 
for a sufficient number of data points for the regression analysis, the 
matchup restrictions are somewhat relaxed from those in the Northern 
Hemisphere. For each 24-hour period surrounding the 0000 GMT and 1200 
GMT radiosonde times, all radiance data is placed into a five degree 
latitude by five degree longitude grid for the latitudes 20®S to 80°S. 
The radiance fields are analyzed and the radiances in each of the eight 
channels are linearly interpolated from the gridded data to the station 
locations. The data match-up is restricted to times and locations when 
there is original satellite data within a geographic box 600 nautical 
miles on a side surrounding the station. This technique is used for 
all levels. 

The -election of the regression equations to be used in computing 
thicknesses is based on certain statistical and practical criteria. The 
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step-wise, multiple regression technique uses an F test to determine if 
additional independent variables should be introduced into the particular 
regression under consideration. In addition, the number of independent 
variables is restricted to a maximum of four. Most of the regression 
equations have three independent variables; however, a few have less 
than three. The standard error of estimate for each of the eight layers 
is given in Table 1 for each hemisphere. These statistics are for the 
three latitude bands combined (area weighted), and are for the dependent 
data on which the regression equations are based. 

No regression equations could be determined for the 10-30 mb layer 
in the Southern Hemisphere because of a complete lack of data above 
30 mb. In the analysis that follows the Northern Hemisphere 10-30 mb 
regression equations are applied to the Southern Hemisphere radiances 
to produce thickness fields in the 10-30 mb layer in the Southern 
Hemisphere. 
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Table 1. Hemispheric standard errors of estimate for thickness 


regression equations 
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Geopotential heights for the SIRS-based data set are determined by 
summing the satellite-based thicknesses from a conventional 1000 mb height 
field. In the Northern Hemisphere the NMC 1200 GMT 1000 mb height fields 
are used. In the ^ uthem Hemisphere 1000 mb height fields are obtained 
by converting daily surface pressure analyses of the Commonwealth Bureau 
of Meteorology, Australia. Tlie Southern Hemisphere surface pressure 
analyses are nearly hemispheric in scope, but where blank they are 
filled in subjectively by a combination of extrapolation and climatology. 
Geopotential heights are therefore determined at the following levels: 

700 mb, 500 mb, 300 mb, 200 mb, 100 mb, 50 mb, 30 mb, and 10 mb. In 
this study only thickness (temperature) information above 100 mb and 
geopotential height information at and above 100 mb is used. 


3. Technique used to calculate vertical motions 

The technique used here to calculate the zonal ly-averaged vertical 
motion io similar to that used by Mahlman (1966) for the Northern 
Hemisphere with conventional data. It is basically the thermodynamic 
technique integrated around latitude circles. The total derivitive of 
• tential temperature (0) ’s expanded to give 


de 90 
3t = Tt 


90 
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where is vector horizontal wind on a surface of constant pressure, and 
0 ) is the vertical velocity in pressure coordinates. Expansion of the 
advection term in equation 1 gives 
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Rearrangement of the terms gives 


r 3PJ 
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at 


dt 
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(3) 


Next, equation 3 is integrated over the area bounded by two latitudes 
4>x and ij> 2 , in this study usually lO” apart. 



+ i |cos<^2 [v^^e] - cos<^l } 

X X 

’^2 

[ev.Vp]^ cos(j) d^ , (4) 

where a is the Earth's radius and [ ]^ implies averaging around a 
latitude circle. 

The averaged quantities in the last two terms of equation 4 are 
expanded to give 
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where 

(7) 


Substituting equations 5 and 6 into equation 4, and using 7 results in 
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In equation 8 the second and third lines represent the effect of the 
divergence of the eddy flux and of the mean motion flux, respectively. 
The fourth and fifth lines represents the effect of divergence. 

The last term in equation 8 involves the correlation of 9 and V*vp 
and is assumed small. The second-to-last term in equation 8 is well 
approximated by 
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where [e]^ ^ is the average potential temperature for the area between 
4ii and <^2* The right side of equation 9 is integrated to obtain 
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Neglecting the last term in equation 8 and using the right side of 
equation 10 to replace the second-to-.'ast term gives 
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The last two lines in equation 11 are of opposite sign and are nearly, 
but not quite, equal. With a temperature gradient of 10®K/10® latitude 
and with 9 approximately 20C®K, [9]^ is different fror either 9 at 
^2 or 6 at (^2 only about two or three percent. Thus the last two 
lines in equation will tend to cancel and their sum is negligible. 

The integrals in equation 11 are evaluated, for example, by the 
approximation 


[ll] COS(J. d.J,= cos (<^2 _ (12) 


where 



*1 + < t >2 

is evaluated at latitude s— — ■ 


With the last two lines 


eliminated as negligible and with the approximations exampled by equation 


12, equation 11 becomes 
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Equation 13 is -ised to calculate the zmially-averaged vertical action. 

4. Scrtherr. and South err. Hemisphere Vertical aotion patterns 

The static stability, local tcaperature chcnge and heat flux conver- 
gence can be evaluated froa the satellite-derived thickness anc height 
fields. The local change tera is of the order of 0.1*iC/day. while the 
heating due to the flux convergence is of the order of l*K/day in polar 

liititudcs. The diabetic term is evaluated by using heating rates 

dl6 

conputed by Rodgers (1967) froa cliaatological data. The vera ^ is 
negative everywhere duiing the polar night in the stratosfdierex and is 
cf the order of l"K/day. Because the diabati; heating tera varies less 
with latitude than the flux divergence tera, the latter is aost iiq>ortant 
in determining la*'tudinal variability of u. For this reason the 
enthalpy flux as a function of latitude is ex-^ined. 

Figur.' 1 shows the poleward heat flux in the 10-30 ah layer averaged 
over January 1970 for the Northern Hatisidiere and over July 1969 for the 
Southern Hemisphere. Both SIRS and NMC curves are presented for the 
Northern Hemisphere. The geostrophic approximation is used for these 
computations. In both hemispheres the peak poleward flux is reached at 
60®. It should be remesibered that the polar vortex in the Northern 
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Fig. 1. Poleward heat flux in the 10-30 nb layer-Jaimary 1970 in the 
N. H., July 1969 in the S. H. 
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Fig. 2. Poleward heat fluxes in the stratosphere of the S, H. for 
July 1969. 
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Heaisphere broke down in early January that year. Thus the Northern 
Hesis^diere curves aay not be represottative of normal (non-breakdown) 
conditions. The two SIRS curves peak at alaost exactly the same value. 
The total flux convergence between 60* and the pole is about equal in 
the two hemispheres. However, the shape of the two curves is different. 
Hie Southern Hemisphere curve shows the flux decreasing very rapidly just 
poleward of 60*S and gradually levelling out. The greatest flux con- 
vergence is in the area between 60*S and 75*S. The Northern Hemisjdiere 
curve has a relatively constant decrease with latitude, indicating that 
the flux convergence ^ > spread out more uniformly than that in the 
Southern Hemisfdiere. This difference in shape of the enthalpy flux 
curves is also evident in the 30-50 mb and 50-100 mb layers, but it is 
not so obvious. 

The enthalpy flux distributions for all three layers for the 
Southern Hemisphere July 1969 are given in figure 2. In all three 
layers there is a rapid decrease in the flux poleward of the peak, then 
a flattening out of the curve in very high latitudes. The poleward 
enthalpy flux at 80*S in the 50-100 mb layer is slightly negative, 
meaning that between 80*S and the pole there is a flux divergence. The 
flux in the 30-50 mb layer at 80*S is very small positive, so that in 
this layer the flux divergence term in equation 13 is very small. In 
the Northern Hemisphere there is substantial flux convergence between 
80*N and the pole (see figure 1) . Thus there is a significant differ- 
ence between the hemispheres in the patterns of heat flux and heat flux 
convergence in the polar-night stratosphere. 

The heat fluxes, the local change in temperature and the static 
stability can be determined from the SIRS-based structure. The diabatic 
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heating, is essentially entirely due to radiative heating or cooling, 
and can not be calculated without a knowledge of the vertical distribu- 
tion of ozone. Even if the ozone distribution is known, the calculation 
is difficult. Rodgers (1967) has calculated tropospheric and strato- 
spheric radiative cooling rates fron cliaatological data in the Northern 
ttoaisi^ere. No such c<Mputations are available for the Southern Hemi- 
sjdiere. Table 2 gives the values of ^ (adapted froa Rodgers (1967)) 
used in this study as a function of latitude and layer. Identical values 
are used in each heaisphere. Rodgers only extended his calculations 
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Latitude 
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O 

00 

O 

O 

60* 

50* 

50* 

10-30 mb 

-1.2 

-1.2 

-1.0 

-0.8 

-0.5 

30-50 mb 

-0.9 

-0.9 

-0.8 

-0.6 

-0.3 

50-100 mb 

-0.6 

-0.6 

-0.5 

-0.4 

-0.2 


Table 2 Values of ^ used in vertical aoticm calculations. 

(Based on Rodgers (1967) calculations) Units: *C/day. 

to 70"N. This author assumes that the values at 70* hold at all 
latitudes poleward of 70*. The values in Table 2 for the 10- SO ad> 
layer are frcs Rodgors* calculations at 20 gd>, the values for the 
30-50 mb layer are an average of the 30 mb and 50 mb calculations, and 
the values for the 50-100 mb layer are from Rodgers' calculation for 
70 mb. In applying these cooling rates to the Southern Hemisphere the 
author is assuming that there is no substantial difference in the actual 
rates between the hemispheres. 
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Fhe results of evaluating the terns in equation 13 und solving for 
u in the Northern Henisphere for January 1970 ifith both NMC and SIRS- 
based structure are shoim in figures 3a and b. The raC-base-:^ calculaticn 
C>a) shows losing notion in high latittides and sinking in low latitudes. 
Tlie axis of rising notion slopes poleward with height aiid the naxinun 
'ising notion is located at 20 mb at the pole. The pattern of rising 
ration over the pole and sinking in nid-latitudes is typical of patterns 
found by other authors (e.g., Mahlnan, 1966) for the Northern Hatisphere. 
Mahlnan (1966) also showed that this pattern of rising and sinking notion 
occurs before, during and after breakdowns. Miller, Brown and Canpana 
(1972) using the quasi-geostrophic, adiabatic onega equation obtain a 
cross-section of u for the period 1-15 January 1970 «diich also Indicates 
the sane general pattern of ascent and descent as in figure 3a. Thus 
the use of equation 13 and NMC structure appears to produce a cross- 
section of vertical notion in the stratosphere which closely resenbles 
previous estimates made for different years, and also resenbles a 
calculation made for the sarc year. 

Figure 3b is the vertical notion pattern for the Northern Henisidtere 
calculated ^.sing the heat fluxes, local te^>erature changes and static 
.‘:tability derived from the SIRS-based structure. There are some simi- 
larities and some differences between figures 3a and b. In figure 3b 
there is rising motion over the pole in the higher layers and sinking 
motion in rldlsticudes. There is also a small area of rising motion at 
60*1! at mb. The maximum rising motion is found at or near the pole 
and in the top layer (10-30 mb). However the volume of rising notion 
is much less than that shown in figure 3a. Magnitudes of maximn 
rising and sinking motions are underestimated. 
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Fig. 3. Zonal ly-averaged vertical Mtion patterns for the N. H. for 
January 1970. Units: mb/day. Stippled area denotes area 

of ascent. 
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The main factor in the difference between figures 3a and b is the 
underestimation of eddy heat fluxes in the SIRS-based calculation. This 
underestimation is evident in figure 1. Although the SIRS-based curve 
underestimates the magnitude of the actual fluxes « it reproduces the 
shape of the NHC curve and the locations of flux divergence and conver- 
gence. This is also true for the two lower layers (30-50 mb and SO-100 
mb). The underestimation by the SIRS-based fluxes varies with latitude 
and between layers. The minimum underestimation between 45*N and 7S*N 
is 56 per cent. To at least partially correct for the underestimation^ 
a factor of 1.8 (based arbitrarily on the 56 per cent figure) is applied 
to the SIRS-based flux values and the vertical motion c(Mq>utation redone. 
The results are shown in figure 3c. The general rising motion in high 
latitudes at high levels and the sinking motion in mid-latitudes are 
still present. The region of rising motion has expanded and indicates 
the slope with height of the lAIC-based cross-sect im. Magnitudes of 
maximum rising and sinking motion are also in closer agrentent. The 
adjustment factor of 1.8 as based on the minimum underestimation is 
chosen because the author believes that the larger underestimations 
present, especially in the lower layers, are due to the low density of 
satellite radiance data in the Northern Hraisphere for January 1970 
relative to the density of observations over the Southern Hwaisphera 
in July 1969, the Southern Hemisphere winter month to be examined. 

The results of the vertical motion calculations for the Southern 
Hemisphere for July 1969 are shown in figures 4a and b. Part a of the 
figure is based on fluxes calculated directly from the SIRS-based 
structure and part b is based on calculations n:?de with the fluxes 
adjusted by the factor of 1.8. Both parts of figure 4 show an area of 
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Pig. 4. Zonal ly-averaged vertical notion patterns for the S. H. for 
July 1969. Units: mb/day. 
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descending motion centered at 50"S, an area of ascent centered at about 
73**S and area of descent over the South Pole. With the assumption of 
ascent in low latitudes (<30**S), figure 4 indicates a three-celled 
meridional circulation compared to the two-celled structure in the 
Northern Hemisphere. The predominant sinking motion over the Southern 
Hemisphere pole has been previously postulated by Reiter (1969, 1971) 
on the basis of ozone measurements over Antarctica. Since Rubin and 
Weyant (1963) have already noted a general subsidence below 75 nd) to the 
surface over the South Pole it is evident that ozone-rich air in the 
lower stratosphere is transported downward to the surface to produce the 
midwinter peaks in surface ozone observed over Antarctica (Wisse and 
Meerburg, 1969). 

5. Southern Hemisphere variations during July 1969 

The early part of July was characterized by a midwinter minor 
warming in the Southern Hemisphere stratosphere (Fritz and Soules, 1970; 
Miller, Finger and Gelman, 1970). Except for the last few days, the 
remainder of the month was relatively quiet. The vertical oration patterns 
during the minor warming and during the quiet period exhibit significant 
differences. The meridional circulation for 1-10 July 1969 (during the 
minor warming) is given in figwes Sa and b. The pattern for 18-27 July 
1969 (quiet period) is given in figures 6a and b. During 1-10 July the 
descending branch over the South Pole disappears and the ascending and 
descending branches at about 70”S and 50*S respectively have larger 
magnitudes than the monthly averages in those locations. Even though 
the descending branch over the pole disappears, the ascending branch is 
still centered between 70"S and 75"S. During the quiet period (18-27 
July 1969) the vertical motion pattern shows very strong descent over 
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l-»0 JULY 1969 S.H. 
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Fig. 5. Zonally-averaged vertical notion patterns for the S. H. for 
the period 1-10 July 1969, during a ninor nid-winter warning. 
Units: wb/iKy. 
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Fig. 6. Zonal ly-averaged vertical motion patterns for the S. H. for 
the period 18-27 Jul)r 1969, a quiet period. Units: mb/day. 
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the pole (figure 6a and b). According to figure 6b the descending motion 
covered the area poleward of about 78°S with maximum downward velocity 
occurring in the lowest layer of calculation. The ascending branch at 
70**S and the descending branch farther equatorward are less intense than 
the monthly average (figure 4) and much less intense than during the 
minor warming (figure 5). Thus during minor warmings in the Sou' em 
Hemisphere it appears that the descending branch over tne pole u..sappears. 
This disappearance may have consequences in the distribution of total 
and surface ozone over Antarctica. The 18>27 July quiet period is 
probably more typical of normal conditions during the Southern Hemisphere 
winter. 

The difference in vertical motion patterns between the periods 1-10 
July and 18-27 July 1969 is closely related to the eddy heat flux distri- 
butions for those two periods. These distributions are shown in figures 
7a and b. The peak fluxes during 1-10 July are approximately three times 
as great as those for the 18-27 July period. Thus the flux divergences 
on either side of the peak are also about three times as great. In very 
high latitudes another important difference is that although both the 
50-100 mb and 30-50 nd) layers exhibit eddy flux divergence between 80*S 
and 90”S during the period 18-27 July, all layers exhibit flux convergence 
in that location for the 1-10 July period. This is the reason for the 
existence of the difference in vertical motion patterns over the pole 
between the two periods. 

Fritz and Soules (1970) examine a number of Southern Hemisphere 
stratospheric temperature fluctuations during the fall, winter and spring 
of 1969 using SIRS channel 8 data. Between 25 June and 10 July they 
indicate a large increase in channel P radiance in the southern polar 
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regions. Acc{»q)anying the warming (increasing radiance) in the high 
latitudes of the Southern Hemisphere, they also find decreasing radiance 
(cooling) in Southern Hemisphere low latitudes and in the Northern 
Hemisphere. The early July case is just one of several sir ir coiQ)en- 
sating changes noted by Fritz and Soules. In this sectic.i c has been 
noted that there are significant differences in the edd> fl'ix of enthalpy 
patterns and in the mean meridi<mal motion patterns between the time of 
the early July warming (1-10 July) and the quiet period (18-27 July). 
These observations support the Fritz and Soules postulate that irregular 
(non-seasonal) fluctuations in stratospheric temperature are related to 
changes in the large-scale eddies and/or changes in the mean meridional 


circulations 
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6. Conclusions 

In this stiidy zonal ly-averaged vertical actions in the stratosphere 
of both hemispheres during winter hav«> been calculated using atmospheric 
struct e acrived from satellite radiance measurements. The calculations 
reconfirm the work by previous authors that the Northern H^isphere has 
a two-celled structure with rising motion over the pole and at low lati- 
tudes and sinking motion in mid- latitudes. The Southern Kemsiphere 
stratosphere, however, has a three-celled structure with sinking over the 
pole and in mid-latitudes and ascending branches at about at ~0°S and 
in low latittjdes. 

The Southern Hemisphere vertical motion patterns are significantly 
different during a time of a minor midwinter warming and during a quiet 
period. During the minor midwinter warming the descending branch over 
the pole disappears and the ascending and descending branches at about 
70“S and 50“S, respectively are strong. During the quiet period the 
descending branch over the pole is strong while the lower latitude 
branches are relatively weak. 
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